arXiv:1503.04154v2 [cond-mat.mes-hall] 18Jun2015 


Symmetries of quantum transport with Rashba spin-orbit: Graphene spintronics 


Leonor Chico,^ A. Latge,^ and Luis Brey^ 

^Instituto de Ciencia de Materiales de Madrid, Consejo Superior de Investigaciones Cientificas (ICMM-CSIC), 

C/ Sor Juana Ines de la Cruz 3, 28049 Madrid, Spain 
‘^Instituto de Fisiea, Universidade Federal Fluminense, Niteroi, Av. Litordnea sn 24210-340, RJ-Brazil 

(Dated: June 19, 2015) 

The lack of some spatial symmetries in planar devices with Rashba spin-orbit interaction opens 
the possibility of producing spin polarized electrical currents in absence of external magnetic field 
or magnetic impurities. We study how the direction of the spin polarization of the current is related 
to spatial symmetries of the device. As an example of these relations we study numerically the 
spin-resolved current in graphene nanoribbons. Different configurations are explored and analyzed 
to demonstrate that graphene nanoflakes may be used as excellent spintronic devices in an all¬ 
electrical setup. 


I. INTRODUCTION 

One of the most important challenges in physics and 
materials science is the exploration of novel systems and 
physical mechanisms for spintronics, with the aim of de¬ 
signing high-speed and low-power devices.In particu¬ 
lar, the production and detection of spin-polarized cur¬ 
rents by electrical means is a newly explored route to¬ 
wards this goal. Remarkably, spin-dependent transport 
can be achieved in systems with spin-orbit interaction 
(SOI) in the absence of ferromagnetic contacts or exter¬ 
nal magnetic fields.^ The extrinsic Rashba SOI couples 
the orbital motion and the spin of the electron under 
an external electric field, allowing for the manipulation 
of spins without breaking time reversal symmetry.^ The 
possibility to tune SOI in Rashba systems is increasing 
the exploration of materials and devices which exploit 
this effect. 

Devices and materials with Rashba SOI are intensely 
studied for spintronic applications, such as transition 
metal dichalcogenides, (TMC)^’^ surfaces or novel two- 
dimensional (2D) nanosheetand one-dimensional 
systems, i.e., nanowires.Theoretically, spin trans¬ 
port has been studied in semiconductor quantum wires 
with Rashba SOI and magnetic field modulations, 
and there are theoretical proposals for spin filters based 
in 2D graphene^^ and graphene nanoribbons with specific 
geometries.Indeed, the experimental ability to fabri¬ 
cate precise graphene nanoribbons (GNRs) by bottom-up 
fabrication processes^’^^ or by epitaxial growth on silicon 
carbide^^ signals these ribbons as potentially fundamen¬ 
tal building blocks in nanoelectronics and is undoubtedly 
one of the reasons to suggest them as spin-orbit-based de¬ 
vices. 

In pristine graphene the intrinsic SOI is negligible; 
however, other spin-orbit couplings induced by differ¬ 
ent mechanisms such as hydrogenation, chemical func¬ 
tionalization or proximity effect with materials with 
strong SOI, have been theoretically proposed^^^^^ and 
experimentally realized.The experimentally re¬ 
ported enhancement of SOI in graphene due to weak 
hydrogenation,^^ gold hybridization,^^ or proximity to 


WS 2 ^^ is of three orders of magnitude or larger, indi¬ 
cating the possible use of graphene in spintronic devices. 
Recently, a giant spin-Hall effect has been experimentally 
measured in graphene, due to the dramatic increase of 
SOI produced by Cu atoms on graphene grown by chem¬ 
ical vapor deposition, reporting SOI splittings around 20 
meV.^^ Intercalation of An atoms in graphene grown on 
Ni has lead to SOI splittings around 100 meV due to 
hybridization with gold atoms.Calleja et al. report 
larger values when Pb is intercalated between graphene 
and the Ir substrate.In addition, spin angle-resolved 
photoemission spectroscopy experiments suggest that a 
large Rashba-type SOI can be tuned in graphene by the 
application of an external electric fieldin samples 
grown on Ni(lll), splittings larger than 200 meV have 
been reported.In nanotubes and curved graphene, hy¬ 
bridization between tt and a bands induces a SOI effect 
larger than in graphene,so folds and wrinkles may also 
increase SOI in graphene systems.^^ Therefore, altering 
graphene by hybridization is an active route to achieve 
SOI values of interest for spintronics. 

In this work we show that the existenee or ahsenee of 
spin-polarized eurrents in planar deviees ean be predieted 
based on a eombination of time reversal and spatial sym¬ 
metries. We give a global and comprehensive vision of 
the possible spin polarization directions for planar ge¬ 
ometries in Rashba systems. As an interesting example, 
we discuss here the spin-dependent conductance of dif¬ 
ferent GNRs, which also present electron-hole symmetry. 
We show that anti-zigzag and anti-armchair GNRs with 
an electric-field-induced Rashba coupling in the central 
region, can produce a spin-polarized current in the di¬ 
rection perpendicular to the ribbon, without the need of 
breaking time reversal symmetry. Furthermore, we also 
demonstrate that spin-polarized currents in the transver¬ 
sal direction can appear in any nanoribbon of constant 
width. 

Indisputably, the discussed effects will be larger in 
other materials, such as TMG^ or germanene and 
stannene nanoribbons,^ but we choose graphene as a pro¬ 
totype system to study SOI effects, which can be modeled 
with a simple Hamiltonian,^^ with the idea that sym- 
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metry reasoning is equally applicable to other materials 
which require a more sophisticated calculation. Notwith¬ 
standing, given the giant SOI experimentally reported 
in graphene,we propose the use of perfect graphene 
nanoribbons as magnetic-field-free spintronic devices. 


II. GEOMETRY AND DEFINITIONS 


We consider planar two-terminal devices as the one 
shown in Fig. 1. The current flows between left (L) 
and right {R) contacts along the longitudinal -direction. 
The spin-orbit coupling occurs in the central part of the 
device that is perturbed by a Rashba-like SOI generated 
by an electric field applied perpendicular to the ribbon; 
i.e., in the ep-direction. The unitary vector Ct defines the 
transverse direction of the ribbon, i.e., across its width. 
The conductance {E) indicates the probability that 
an electron in the left electrode with energy E and spin 
pointing in the cr-direction reaches the right electrode 
with spin pointing in the a'-direction. 



FIG. 1. (Color online) Schematic drawing of the device geom¬ 
etry. Left (L) and right (R) contacts are pristine nanoribbons 
without SO interaction. The conductor (C), shaded in red, is 
the central part of the device with Rashba SOI, length i, and 
width W. 


III. SYMMETRY CONSIDERATIONS 


In planar quasi-one-dimensional devices as those con¬ 
sidered here, one can expect the following spatial sym¬ 
metries: 

(i) C 2 rotation around Cp. Under (^ 2 , the spatial, mo¬ 

mentum and spin components change as {ei^Ct^Cp 
-ei, -et,ep), {pi,Pt,Pp -> -pi, -Pt,Pp), and {ai,at,ap 
—(Ju cTp). Therefore, the conductance of the device is 
invariant under these operations. This amounts to inter¬ 
change the left and right electrodes, and invert the spin 
direction along the ei or Ct directions, i.e., = G§^,, 

where a indicates a spin projection opposite to a, and 
cr'=±cr. For the spin direction perpendicular to the de- 
vice, we get 

(ii) Longitudinal mirror symmetry M/. For M/, 
the spatial and momentum components transform as 


{ei,et,ep ei,-et,ep) and {pi,Pt,Pp ^ Pi,-Pt,Pp), re- 
spectively. For the spin components, recalling that the 
spin transforms as an axial vector, we have (a/, (Jt^cFp 

—CTp)- Thus, this symmetry does not change the 
roles of the electrodes, but changes the sign of the spin 
projection when the spin direction is either ei or e^, lead¬ 
ing to the relation = G^^,. Notice that this sym¬ 

metry does not give any relationship for the conductance 
when the spin direction is along Cf. 

(hi) Transversal mirror symmetry Mf. Under 
the spatial and momentum components transform as 
{euet.Cp -euet.Cp), {puPt.Pp -p/,pt,Pp), whereas 
the spin changes as {(Ti^cjt^cFp —at^ —cTp). There¬ 

fore, the relation G^^, = G^^, is obtained when the spin 
is pointing in the Ct or Cp directions. Otherwise, for the 
spin pointing in the longitudinal direction e/, we obtain 
the relation G^^, = G ^^,. 


TABLE 1. Symmetries and the corresponding conductance 
relations derived from them. First column: symmetries; sec¬ 
ond column: spin projection directions; third column: spin- 
resolved conductance relations. In this Table a'^Ea. 


0 

Cl, Ct, Cp 

/^LR /^RL 

^acr' — 

U 

Clf Ctf Cp 


Ml 

Cl , Cp 

/^LR /^LR 

Mt 

ei 

^LR ^RL 

Ct 5 Cp 

/^LR /^RL 

C 2 

ei,et 

^LR ^RL 

Cp 

^LR ^RL 


Besides spatial symmetries, in the absence of magnetic 
field time reversal symmetry (0) holds, implying G^^, = 
G§/^. Finally, in presence of electron-hole symmetry (U), 
the conductance as a function of the energy E satisfies 
G^^,{E) = G^f^{—E). All these relations are gathered 
in Table 1. These symmetry rules allow us to predict 
the possibility of obtaining spin-polarized currents in any 
planar devices. For an incident unpolarized current from 
the left electrode, the spin polarization of the current in 
the right electrode in the 5-direction (5 = /, t,p) is defined 
as 


P I /^LR I /^LR /^LR /^LR /-|\ 

s ~ ^ss T ^ss ^SS ^SS • GJ 

Therefore, from symmetry considerations we obtain that 
in systems with Mi symmetry the spin polarization of the 
current in the ei and Cp spin-directions are zero, whereas 
in all the other cases a spin polarized current can be 
obtained, albeit with different intensities and for various 
reasons. 
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FIG. 2. Schematic drawing of the four considered central 
conductors and their relevant spatial symmetries. 


IV. SPIN-POLARIZED CURRENT IN 
GRAPHENE NANORIBBONS 

We consider the simplest possible geometry: a 
graphene nanoribbon where the central region has a finite 
Rashba SOI. The length of the conductor i is given by 
3accN for an armchair (AC) GNR and by \/3accN for a 
zigzag (ZZ) one, with N being the number of longitudinal 
unit cells and Ucc the carbon bond length in graphene. 
The ribbon width, VF, is defined by the number of dimers 
(zigzag lines) across the width of the armchair (zigzag) 
nanoribbon, given hereafter by M. As the focus is on 
transport properties, in the case of AC terminations we 
restrict the study to metallic ribbons, i.e., those with 
M = 3q 2, q being an integer. 

Infinite zigzag GNRs have transversal mirror symme¬ 
try Mt for any M, but this symmetry is not present in 
a rectangularly cut finite-size fiake of length V, such as 
those depicted in the upper panels of Fig.2. The parity 
of M defines two kinds of zigzag flakes: even-M zigzag 
GNRs that have a longitudinal mirror symmetry Mi and 
odd-M zigzag GNRs which are invariant under (72- The 
M-odd zigzag ribbons are commonly called anti-zigzag 
(AZZ), as shown in Fig. 2. For armchair GNRs, the more 
symmetric configurations, with M/, Mt and C 2 symme¬ 
tries, happen for M odd. M-even AC GNRs only have 
Mt, both in the infinite case and in the finite flake and 
they are called anti-armchair (AAC) ribbons, see Fig. 2. 

The application of the symmetry relations summa¬ 
rized in Table I to the GNR flakes depicted in Fig. 2 
tells US the possibility of obtaining spin-polarized cur¬ 
rents in graphene-based devices. Fig. 3 shows the rela¬ 
tion between the spin-resolved conductances for different 
graphene flakes. Note that for Mt and C 2 symmetries, 
it is necessary to consider also time reversal symmetry 
0 to obtain the relations between the spin-polarized cur¬ 
rents flowing in the same direction, as indicated in Fig. 3. 
The same symmetry arguments can be easily generalized 
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FIG. 3. Graphical summary of our results for graphene 
nanoribbon flakes. For each case, the first column lists the 
spin direction; the second column shows the corresponding 
spin-resolved conductance relation with the spatial symme¬ 
try employed for its derivation, indicating also whether time- 
reversal symmetry is needed (with a clock icon) in order to 
obtain the final result shown in the third column, i.e., which 
spin conductances are equal or not. Up and down arrows are 
referred to the spin projection direction indicated in the first 
column. In the last column the red (blue) color of the rectan¬ 
gle surrounding the final spin conductance relations indicates 
the possibility (impossibility) of getting spin-polarized trans¬ 
port. 

along the lines described in this work to other graphene 
flake geometries. However, the use of the symmetry re¬ 
lations does not give ns information on the magnitude of 
the spin polarization. In order to quantify the spin polar¬ 
ization it is necessary to perform microscopic calculations 
taking into account all the details of the discrete lattice 
and distinguishing between different symmetries. To do 
that we use a tight-binding formalism for computing the 
conductance of different GNRs. 


A. Microscopic calculations 

The Rashba spin-orbit interaction in graphene can be 
described in the nearest-neighbor hopping tight-binding 
approximation^^’^^. The total Hamiltonian is M = 
Hq + where Hq is the kinetic energy term, Hq = 
with t being the nearest-neighbor hopping 
and c[(x destruction and creation operators for 
an electron with spin projection a in site i and j, respec¬ 
tively. The Rashba SOI contribution is given by 

4a [(o- X dij) ■ ep]al3Cji3 , (2) 

“ <i,j> 

a,/3 

with Xr being the Rashba SOI strength that can be tuned 
by the electric field intensity, a are the Pauli spin matri¬ 
ces, dij the position vector between sites i and j, (3 
are the spin projection indices. Due to the analytical ex¬ 
pression for Hr (Eq. 2), a new symmetry occurs. If Xr 
changes sign, the spin polarizations on the longitudinal 
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and transversal directions also change sign, whereas the 
polarization in the perpendicular direction is unaffected, 
i.e., P,(-A) = -Ps{X),{s = = P,{X), {s = 

p). Notice that the sign of Xr is determined by the sense 
of the electric field which originates it. 

We consider a graphene device composed of a cen¬ 
tral fiake with Rashba SOI and two semiinfinite pris¬ 
tine nanoribbons of the same width as the conductor, 
see Fig. 1. The conductance is computed in the Kubo 
approach by using the Green function formalism^^’^^. 
The spin-resolved conductance is given by = 

, where is the advanced (re¬ 
tarded) Green function of the conductor and T^^^^ = 
ilYl^LiR) o- ~ "^L{R) cr] written in terms of the L (R) lead 
selfenergies a- equivalent expression can 

be reached from a scattering approach, equivalent to the 
Kubo formalism, which evidences more clearly that the 
symmetry of the conductance matrix is the same as that 
of the Hamiltonian of the system. 


B. Numerical results 

In order to verify the previous symmetry analysis, 
we perform numerical calculations for several graphene 
nanoribbons of similar widths and lengths, but for dif¬ 
ferent symmetries and spin polarization directions. The 
most common expression for the Rashba Hamiltonian 
chooses z as the electric field direction, and therefore sets 
this as the spin quantization axis, even though it is not 
the most favorable from both the quantitative and the 
symmetry viewpoint. For some systems, such as TMG 
(WSe 2 and M 0 S 2 ), h has been experimentally reported 
that spin polarization can occur in this particular setup. 
Indeed, it has been dubbed a ” Zeeman-type” spin split¬ 
ting with an electric field, due to the fact that the spin 
and the external field are in the same direction.^ We first 
analyze the results for this spin polarization direction, 
and then concentrate on the best configuration for the 
obtention of a spin-polarized current, namely, with the 
spin transversal to the current and the electric field. 


1. Spin direction perpendicular to the plane of the ribbon. 

When the spin direction is along e^, fiakes AG and 
ZZ which have Mi symmetry cannot present a net spin- 
polarized current. We have checked numerically this re¬ 
sult, not shown here. On the contrary, the anti-armchair 
and anti-zigzag GNRs do not show Mi symmetry, so they 
may have a net spin-polarized current. Indeed, we have 
found in both cases a finite spin-polarized transport; how¬ 
ever, in this geometry the effect is small, especially for the 
anti-zigzag fiakes. Therefore, for the sake of clarity of the 
spin-dependent conductance plots, we use a large value 
of SOI for this configuration, = 0.3t. Fig. 4 shows 
the spin-resolved conductances for an AAG graphene 
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FIG. 4. Spin-resolved conductances and spin polarization of 
the current as a function of the energy for an 8-AAC GNR of 
length N = 6 for spin direction perpendicular to the plane of 
the GNR and Xr = 0.3t. Panel (a) shows the spin-conserved 
and panel (b) the spin-flip conductances. Panel (c): Resulting 
spin polarization of the current as a function of the energy. In 
(d) we present the corresponding band structure of the infinite 
ribbon with Rashba SOI (thick blue lines) and without SOI 
(thin black lines). The top left schematic drawing indicates 
the spin projection direction used for this plot, shown with a 
thicker arrow. 


nanoribbon with M = 8 and length N = 6. The AAG 
fiakes have Mt symmetry and therefore = G^^^ 

so a net spin-polarized current occurs because the spin- 
conserved conductances are different, G^^ 7 ^ 

In contrast, as the AZZ case has C 2 symmetry, the 
spin-conserved conductances are equal and the spin po¬ 
larization occurs because Fig. 5 presents 

the spin-resolved conductances for an AZZ ribbon of 
width M = 9 and length N = b. The spin current is 
smaller in this case. We attribute this to the fact that 
the infinite AZZ ribbon has also Mt symmetry, absent 
in the fiake, which would yield equal spin-flip conduc¬ 
tances. Therefore, in this case the spin-polarized cur¬ 
rent is clearly a finite-size or scattering effect, due to 
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FIG. 5. Spin-resolved conductances and spin polarization of 
the current as a function of the energy for a 9-AZZ GNR with 
length N = 5 and spin direction perpendicular to the plane 
of the GNR and Xr = 0.3t. (a) Spin-conserved conductances, 
which are equal due to C 2 symmetry; (b) spin-flip conduc¬ 
tances, which give rise to the spin-polarized current, (c) Spin 
polarization of the current, (d) Band structure of the cor¬ 
responding infinite ribbon calculated with Rashba SOI (thick 
blue lines) and without SO interaction. The top left schematic 
drawing indicates the spin projection direction considered for 
this plot, shown with a thicker arrow. 


the boundary between the leads without Rashba and the 
flake with SOI. 

In these two instances we can observe another rela¬ 
tion for the conductances and polarization derived from 
both time reversal and electron-hole symmetries. It can 
be easily verified that the combination of 0 and U (see 
Table I) yields In terms of the 

polarization, this means Ps{E) = —Ps{—E). Figs. 4 (a) 
and 5 (b) are non-trivial examples for the spin-conserved 
and spin-flip conductances, whereas Fig. 4 (c) shows this 
symmetry in terms of the current polarization Pp{E). 

Notice that there is a region around Ep where the spin 
polarization is zero. It corresponds to the energy range 
with only two bands in the energy spectrum, as it can be 


seen in panels (d) of Figs. 4 and 5. In order to have a net 
spin current, more than two bands should be available in 
the system. Using wider ribbons lowers the energies of 
these additional bands. 


2. Spin projection parallel to the transversal direction of 
the nanorihbon. 

Now we discuss the case with spin projection parallel to 
the nanoribbon transversal direction Cf. In this setup the 
spin-polarized conductance is the largest, so this config¬ 
uration is the most relevant from the experimental view¬ 
point. This result can be inferred from the structure of 
the Rashba Hamiltonian. In the continuum model the 
Rashba term takes the form Hr oc {a x k) • E, being E 
the applied electric field. Thus, a maximum contribution 
can be expected when the directions of the electric field 
E^ the current direction k and the spin are orthogonal, 
as it happens when the spin is pointing in the transver¬ 
sal direction. For this spin direction, the symmetries do 
not impose any condition to the spin-conserved conduc¬ 
tances in all considered cases, i.e., G^^ ^ Actually, 

polarized currents are obtained in all four GNRs studied. 

Results for each one of the symmetries illustrated in 
Fig. 2, namely, zigzag (M = 8), armchair (M = 11), 
anti-zigzag (M = 9) and anti-armchair (M = 8) GNRs 
are presented in Fig. 6. The chosen ribbons have similar 
widths and the same fixed length, N = 4. For this spin 
direction, e^, the magnitude of the polarized current is 
similar in the four cases. Moreover, the spin polarization 
of the symmetric GNRs presents magnitudes similar to 
those found in the asymmetric cases. This result may be 
understood on the basis that Mi does not play any role 
for this particular spin orientation. Also, it is evident the 
relation Ps{E) = —Ps{—E), imposed by electron hole- 
symmetry, as given in Table I. 


3. Final remarks. 

Finally, we would like to mention that if the spin pro¬ 
jection is taken along the longitudinal direction, the effect 
is small, although not zero. The particular conductance 
relations are also collected in the graphical summary pre¬ 
sented in Fig. 3. 

As our main focus in this work was the role of sym¬ 
metry in the spin-resolved conductances, we have cho¬ 
sen small flakes, for which the effects are clearer. For 
longer flakes, the most notable (and obvious) difference 
would be that the number of conductance and polariza¬ 
tion oscillations increases, due to the appearance of more 
quasi-localized states. 

It is interesting to mention that the set of symmetries 
discussed here also plays an important role on the size 
dependence of the spin polarization. As discussed previ¬ 
ously, in some systems, the spin polarized currents arises 
because of a finite-size effect, whereas in others it is not 
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FIG. 6. Left panels: Spin-resolved conductances as a function of the energy for an 8-ZZ, an 11-AC, a 9-AZZ and an 8-A AC 
graphene nanoribbons of length A" = 4, enumerated from top to bottom. In all cases the spin is projected in the transversal 
direction of the CNR, indicated with a thicker arrow in the schematic drawing at the top. Right panels show the corresponding 
polarization currents Pt. The results are obtained with Xr = O.lt. 


so. If the polarization is due to a finite size effect, then we 
expect that its maximum value will eventually decrease 
with size, but in any case it should not grow on average. 
On the other hand, if the polarization is due to a lack of 
symmetry present both in the infinite and in the finite 
case, there should be a non-zero polarization for all sizes. 

As an example we have chosen a 9-AZZ ribbon; the cor¬ 
responding flake has only C 2 symmetry, but the infinite 
ribbon has also M^. For spin projection direction along 
Cp, this implies an extra relation in the spin-resolved con¬ 
ductance (see Table I) that yields a zero polarization cur¬ 
rent. However, for spin projection direction along e^, the 
relation is the same as for C 2 , so we expect the polarized 
current to exist for growing size. These size dependences 
are illustrated in Figure 7, where spin-polarized currents 
are shown for the two discussed spin polarization direc¬ 
tions and two flake lengths with Rashba SOI. Besides 
the aforementioned oscillations due to the larger size, it 
is notable the increase of the polarization of the current 
presented in the bottom panel for the longer flake length. 

As a general result, we have shown that for some par¬ 
ticular graphene flakes (Fig. 2), the use of symmetries 
allows us to elucidate which spin-resolved conductances 
are equal and which are different. In the same way, al¬ 
though not presented here, we can also infer that GNRs 
with symmetric chiral edges, as those obtained by open¬ 
ing carbon nanotubes, will behave as AZZ ribbons, be¬ 
cause of their C 2 symmetry. With respect to the size 
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FIG. 7. Spin polarization currents as a function of the energy 
for 9-AZZ GNRs with different Rashba SOI strength, spin 
polarization directions, and ribbon lengths: (a) N=5 and 20, 
Ai?=0.3t and spin projected in the perpendicular direction 
(Pp) and (b) N=4 and 20, Ai?=0.1t and spin projected in the 
transversal direction (Pt)- 

dependence, for wider ribbons the number of bands in¬ 
creases, and the onset of spin polarized currents happens 
for lower energies because of the availability of more than 
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two spin channels. 

V. SUMMARY 

We have studied the symmetries of the spin-resolved 
conductances in planar devices with Rashba SOL The 
combination of spatial mirror reflections and C 2 rotation 
with time-reversal symmetry leads to specific predictions 
with respect to the possibility of obtaining spin-polarized 
currents in such devices. As an example, we compute the 
spin-dependent transport of graphene nanoribbons with 
an applied electric field in finite region. We have shown 
that spin-polarized currents can be achieved if the spin 
polarization is measured in the transversal direction of 
the ribbon for all the ribbon geometries. Furthermore, 
we have analyzed all the basic symmetries and spin direc¬ 
tions, elucidating which configurations can yield a spin- 


polarized current on the basis of symmetry. The inten¬ 
sity and sign of the Rashba spin-orbit coupling may be 
modified by external electric field, opening the possibil¬ 
ity of building an all-electrical spin valve. Our findings 
can be useful for a smart design of spintronic graphene 
devices, being of general application to other materials 
with Rashba SOI. 
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